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I. S. Bcloborodov, 1 ' 2 A. Glatz, 1 and V. M. Vinokur 1 

1 Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA 
2 James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA 
(Dated: February 1, 2008) 

We study electronic transport properties of ferromagnetic nanoparticle arrays and nanodomain 
materials near the Curie temperature in the limit of weak coupling between the grains. We calculate 
the conductivity in the Ohmic and non-Ohmic regimes and estimate the magnetoresistance jump 
in the resistivity at the transition temperature. The results are applicable for many emerging 
materials, including artificially self-assembled nanoparticle arrays and a certain class of manganites, 
where localization effects within the clusters can be neglected. 
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Arrays of ferromagnetic nanoparticles are becoming 
one of the mainstreams of current mesoscopic physics |l|, 
0, S H ■ Not only ferromagnetic granules promise to serve 
as logical units and memory storage elements meeting el- 
evated needs of emerging technologies, but also offer an 
exemplary model system for investigation of disordered 
magnets. At the same time the model of weakly coupled 
nanoscale ferromagnetic grains proved to be useful for 
understanding the transport properties of doped man- 
ganite systems 0] that have intrinsic inhomogeneities. 
Recent studies showed that above the Curie temperature 
these materials possess a nanoscale ferromagnetic cluster 
structure which to a large extend controls transport in 
these systems @, H, & Ell • This defines an urgent quest 
for understanding and quantitative description of elec- 
tronic transport in ferromagnetic nanodomain materials 
based on the model of nanogranular ferromagnets. 

In this paper we investigate electronic transport prop- 
erties of arrays of ferromagnetic grains (llj near the 
ferromagnetic-paramagnetic transition, see Fig. [TJ At 
low temperatures, T < T c s , the sample is in a so called 
superferromagnetic (SFM) state, see Fig. [TJ set up by 
dipole-dipole interactions. Near the macroscopic Curie 
temperature, T c s , thermal fluctuations destroy the macro- 
scopic ferromagnetic order. At intermediate tempera- 
tures T" < T <T->, where Xjf is the Curie temperature 
of a single grain, the system is in a superparamagnetic 
(SPM) state where each grain has its own magnetic mo- 
ment while the global ferromagnetic order is absent. At 
even higher temperatures, T > T§, the ferromagnetic 
state within each grain is destroyed and the complete 
sample is in a paramagnetic state. We consider the model 
of weakly interacting grains in which the sample Curie 
temperature is much smaller than the Curie temperature 
of a single grain, T c s <C T§. 

We first focus on the SPM state, and discuss a 
d— dimensional array (d = 3, 2) of ferromagnetic grains 
taking into account Coulomb interactions between elec- 
trons. Granularity introduces additional energy parame- 
ters apart from the two Curie temperatures, T® and T c s : 
each nanoscale cluster is characterized by (i) the charging 



T < T c s 




T| < T 








Superferromagnet 


Superparamagnet 


Paramagiiet 




FIG. 1: Top: Sketch of a 3d granular system under consider- 
ation showing the different states at different temperatures: 
For T <T°, where T c s is the macroscopic Curie temperature 
of the system, the ferromagnetic grains (superspins) form a 
superferromagnet (SFM); for T<f < T < T c 9 , where T c 9 is the 
Curie temperature for a single grain, the system is in a super- 
paramagnetic (SPM) state; and above T| the system shows 
no magnetic order. In the SPM state the angle between two 
grains (superspins) is denoted by 9. Note, this is an idealized 
picture, see Bottom: Schematic behavior of the resis- 

tivity, p, versus temperature, T, in the different states in the 
absence (B = 0) and presence (B > 0) of a magnetic field B 
aligned with the magnetization of the SFM, cf . [H, Il2l) . 

energy E c = e 2 /(na). where e is the electron charge, k 
the sample dielectric constant, and a the granule size, and 
(ii) the mean energy level spacing 5. The charging energy 
associated with nanoscale ferromagnetic grains can be as 
large as several hundred Kelvins ll| and we require that 
E c /5 ^> 1. The typical sample Curie temperature, T c s , 
of the arrays we consider (and also of doped mangan- 
ites) is in the range (100 - 200) K, @, [13, til; thus the 
temperature interval T c s < T < E c is experimentally ac- 
cessible. To satisfy the last inequality, the size of a single 
ferromagnetic grain, a, should be less than the critical 
size a c = e 2 /(T c s k). The condition E c /6 ^> 1 defines the 
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lower limit for the grain size: a; = [pi/e l v) 1 K D ~ r ) , where 
v = + i*i is the total density of states at the Fermi sur- 
face (DOS) with v«n being the DOS for electrons with 
spin up (down) and D the grain dimensionality [Tij . 

The internal conductance of a metallic grain is taken 
much larger than the inter-grain tunneling conductance, 
which is a standard condition of granularity. The tun- 
neling conductance is the main parameter that controls 
macroscopic transport properties of the samrjle 
consideration of applications to experiments 
we restrict ourselves to the case where the tunneling con- 
ductance is smaller than the quantum conductance [17|. 
In the SPM state the charge degrees are coupled with 
the spin degrees of freedom; to reflect this connection, 
the tunneling conductance can be written in a form 
g t {9) = t °(l + n 2 cos6>), [ij], where g t ° is the tunnel- 
ing conductance in the paramagnetic state [19j |: II = 
(i/| — v^/iy] + fj.) is the polarization factor of a fer- 
romagnetic grain where 9 G [0, w] is the angle between 
two superspins, see Fig. [TJ The tunneling conductance, 
gt{9), achieves its maximum value for parallel spins, 9 = 
(corresponding to the SFM state). In general, the distri- 
bution of angles 9 is determined by some function f(9) 
which depends on temperature (and on external mag- 
netic field): for T < T c s , in the SFM state f(0) is the 
(^-distribution and for high temperatures, T ^> T c s , it 
is constant [an explicit expression for f(9) is discussed 
below Eq. ||3J)]. We denote averages over angles by 
(. . ,) g = fade... f{9) with (l) e = 1. Using this distribu- 
tion we introduce the averaged tunneling conductance: 



5t (m 2 ) S (g t (0)) e = 5t °(l + n 2 m 2 ) 



(1) 



with the normalized "magnetization" m 2 = (cos9) g , e.g. 
to 2 = 1 in the SFM state and m 2 — > for high tem- 
peratures, T ^> T r s . Note, that in general m 2 is not 
the normalized (absolute value of the) magnetization of 
the sample since it only takes into account the angle be- 
tween two neighboring superspins in the plane spanned 
by them (Fig. [TJ). However, close to T£, we can expect 
m 2 w \M.(T)/M S \ 2 , where M s is the saturation value of 
the magnetization of the sample. Below we first discuss 
the Ohmic transport near T* and then summarize the 
main results for the resistivity behavior in the non-Ohmic 
regime. 

Ohmic transport. — To calculate the conductivity 
for weakly coupled grains in the presence of quenched 
disorder, we start with determining the total prob- 
ability for an electron to tunnel through TV grains 
P{9 1 ,...,9 N ) = n£LiA(0i), whcrc Pi(0i) denotes the 
probability for an electron to tunnel through a sin- 
gle grain i with an angle difference 9i of the magnetic 
moment to the previous grain [lil . [2fj| . The prob- 
ability P(9i,...,9n) has to be averaged over all an- 
gles in order to obtain the total tunneling probability 
Ptotai = (P(8i, ■ ■ ■ , 9 N )) 6u ...,e N = nti Pi(™ 2 ) ■ The lat- 
tcr equality follows from the fact that 7-totai factorizes 



into the individually averaged probabilities Vi(m 2 ) = 
(Pi(9i))9i, [HI- The mechanism for electron propaga- 
tion through an array of grains at low temperatures is 
clastic and/or inelastic co-tunneling. The corresponding 
probabilities in the limit of weak coupling between the 
grains are given by Vf{m 2 ) ~ \g t (m 2 )5]/E c and 
Vl n (m 2 ) ~ [g t (m 2 )T 2 ]/E 2 , respectively, \2J$. Assum- 
ing that all probabilities Vi(m 2 ) are approximately the 
same [H[ for each grain, Vi(m 2 ) = V(m 2 ) for all i, and 
expressing them in terms of the localization length £ (m ) , 
defined by V(m 2 ) = exp[— a/^(m 2 )], we obtain [15| 

e el ^ a/ HE c /g t (m 2 )8], C * a/ ln[E 2 c /T 2 g t (m 2 )}. (2) 

Since the characteristic temperature we consider is of 
the order of the Curie temperature, T ~ T c s , the dom- 
inant mechanism for electron propagation is the inelas- 
tic co-tunneling [22| . Following Mott-Efros-Shklovskii's 
theory [23|, [24{, the conductivity can be written as 
a(T, m 2 ) ~ gtijn 2 ) exp[— r/^(m 2 ) — e 2 / (nrT)}, where the 
tunneling conductance g t (m 2 ) is given by Eq. JTJ) and 
r is the hopping distance. The first term in the expo- 
nent accounts for electron tunneling and the second term 
describes thermal activation necessary to overcome the 
Coulomb correlation energy. Optimizing a(T,m 2 ) with 
respect to the hopping length, r = N ■ a, we obtain: 



o~(T, m 2 ) ~ g° t (l + U 2 m 2 ) exp(-^o(m 2 )/T), (3) 

with T (m 2 ) = T [l - (£o/a)ln(l + II 2 m 2 )] being the 
characteristic temperature scale. Here To = T)(™ 2 = 
0) = e 2 /(K^ ), where £o is the inelastic localization length 
given in Eq. (f2J) with the tunneling conductance corre- 
sponding to the paramagnetic state, gtim 2 = 0) = g®. 
The minimal value of the resistivity in the SFM state 
is determined by the minimal value of the energy scale 
T min = T (l-ln2). 

The behavior of the resistivity p(T,m 2 ), inverse of 
the conductivity er(T, m 2 ) in Eq. ([3]), and the mag- 
netization m 2 (T) are shown in Fig. [2] for the follow- 
ing set of parameters: £_o/a = 1, 7o/T c s = 10, and 
II 2 = 0.3. In order to describe typical experimental 
data [l2l ], the normalized distribution function f(9) was 
chosen to ensure a sharp drop in the magnetization at 
T c s , f[a,x] = [a/arctan(l/a)](x 2 +a 2 )~\ with x = 9 /n 
and a(T) = 10(T/T C S -1) for T > T c s and a(T < T") = 
otherwise. The numerical constant in a(T) was taken to 
produce a drop of m 2 in a temperature region AT* with 
AT* /T* ~ 0.1. Note, that the jump in the resistivity 
does not depend on the precise functional expression of 
/(#); but it is sufficient that m 2 decays rapidly in the 
interval AT* < T*. 

For small polarization factors, Ii 2 m 2 -C 1, the ex- 
pression for the energy scale To(m 2 ), can be written 
as T (m 2 ) = T [l - {£, /a)n 2 m 2 }. As a result, we 
obtain for the conductivity of an array of superspins 
a(T,m 2 ) ~ [1 + 7T n 2 m 2 ]exp(- v /T /T), where 7T is 
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FIG. 2: Solid line, right axis: Plot of the normalized resis- 
tivity p(T,m 2 ) / pspm, inverse of the conductivity a(T,m 2 ) in 
Eq. ([3j), vs. temperature for the following set of parameters: 
Co/a = 1, To/T c s = 10, and II 2 = 0.3. Here p S PM = p(Ti,m 2 ) 
with Ti > T* being the temperature at which p(T, m 2 ) is 
maximal and m 2 = m 2 (Ti). Dashed line, left axis: Plot of 
the "magnetization" m 2 versus temperature using an angular 
distribution function defined in the text, 
a temperature dependent function. For temperatures 
T s c < T « T it is given by 7r (Co /2a) ^JT. The 
last expression can be written in terms of the optimal 
hopping length, r opt = £o^T a /T, as 7 T ~ r opt /2a. Since 
the hopping length r opt depends on temperature even 
within the SPM state, where the temperature T satis- 
fies the inequality Xj? < T < E c , one expects to observe 
two different regimes: at temperatures T c s < T < T* 
the dominant mechanism for electron propagation is vari- 
able range hopping [in this regime r opt > a ], while for 
T* < T < E c electrons hop between the nearest neigh- 
bor grains only [r opt ~ a]. The separating temperature 
T* can be estimated using the condition r opt ~ a which 
gives T* ~ T (Co/a) 2 , i.e. T*<T . 

To calculate the magnitude of the jump value of the re- 
sistivity at the transition from the SPM to the SFM state 
in Fig. [2j we introduce the dimensionless resistance ratio 
Ap/p = [pspm - /Osfm]/pspm, where psfm is calculated 
at temperature T = T r s , using the inverse of Eq. ^ , for 
magnetization to 2 = 1; psvm bas to be evaluated at some 
temperature T\ > T c s at which p(T,m 2 ) is maximal and 
magnetization m 2 = m 2 (Ti). The temperature T\ is of 
the order of T c s + AT C S ; but since the magnetization drops 
quickly above T c s , i.e. AT C S <C T*, it is sufficient to take 
Pspm at 7\ ~ for calculating the jump. Using Eq. © 
we obtain the following result: 

Ap ^ 1 + n 2 m 2 jT ( m >)/T c - + y /ThW7Ti U) 

p - i + n 2 ' { ) 

For small polarization factors, II 2 <§C 1, this reduces to 
Ap/p ~ (Co /2a) ^T /T c s IP 2 (1 - to 2 ), and can be ex- 
pressed in terms of r op t as Ap/p ~ IT 2 r op t(T c s )/a. There- 
fore the larger the hopping length, the bigger is the re- 
sistivity jump between the SPM and SFM states. From 
Eq. ((3]) follows that the resistance ratio Ap/p increases 



for small Curie temperatures, Xj? <C Tq. We now esti- 
mate the jump magnitude in Eq. ([4]): Using the realistic 
values: IP 2 = 0.3, m\ = 0.1, Co/a = 1, T /T c s = 10, we 
obtain Ap/p ~ 0.4, corresponding to a jump 7 = 60%, 
where 7 is defined by pspm = (1 + 7) Psfm- This estimate 
agrees with the plot in Fig. O 

Non-Ohmic regime. — So far we discussed the Ohmic 
regime in the absence of an additional external electric 
field (or applied voltage) only. In the presence of an elec- 
tric field E, the hopping conductivity in the paramag- 
netic state is a ~ exp [-r/C - e 2 /{nrT) + eEr/T], [ia j. 
with the inelastic co-tunneling localization length C m = 
a/ln[£ 2 / 5t °(T 2 + (e£a) 2 )], For sufficiently high 

electric fields E > T/eC the tunneling term, exp(— r/C), 
in the conductivity is not important. As a result the 
optimal hopping distance r opt (E) ~ ^^JE^/E, with the 
characteristic electric field E^ = e/(ftC 2 ), and the resis- 
tivity p ~ exp[r opt (i?)/C] are temperature independent. 

Including the magnetization dependent tunneling con- 
ductance g t (m 2 ) in the above consideration, one finds 
that the conductivity in the SPM state in the presence 
of a strong electric field <r(E,m 2 ) is given by Eq. (|3|) 
with %(m 2 ) — > £a(m 2 ) and T — > E. where £o(m 2 ) = 
Eq[1 — (2Co/a)ln(l + II 2 to 2 )] is the characteristic elec- 
tric field, E = To/eCo, and Co = a/ \n[E 2 / ([eaE} 2 g°)}. 
Equation ([3]) with the above substitutions holds for elec- 
tric fields T/[eC(m 2 )] < E < £ {m 2 ). The last in- 
equality means that the optimal hopping length r opt (E) 
is larger than the size of a single grain, a, while the 
first inequality ensures that the electric field E is still 
strong enough to cause non-Ohmic behavior. Using typ- 
ical values |2j for a ps lOnm, To w 10 3 K, n ps 3, and 
temperature T m 10 2 K we estimate the window for 
electric fields as 10 3 V/cm < E < 10 5 V/cm. The re- 
sistance ratio Ap/p in the presence of a strong elec- 
tric field is still given by Eq. (J4| with the substitution 
7o(to 2 ) — ► £o(m 2 ). For small polarizations, U 2 ^ 1, one 
obtains Ap/p ~ (Co/a) ^E /(T^/e^) n 2 (1 - to 2 ) . 

Discussion. — Past experimental studies of self- 
assembled ferromagnetic arrays were dealing either with 
their thermodynamic properties 0, H, 0, 0] , or with do- 
main wall motion (26[. Investigations of the electronic 
transport and magnctoresistance (MR) were mostly re- 
stricted to the SPM state [2?], [IH , where variable range 
hopping was observed. The crossover region near Tj? was 
only studied by numerical methods in the context of man- 
ganite systems @, [|1 [l(| ■ The resistivity dependence be- 
low the Curie temperature T c s presented in Fig.^is differ- 
ent from the schematic behavior shown in Fig.[T]rcflecting 
the experimental data on manganite systems of Ref. [l2| . 
In these materials the resistance in the SFM state (below 
T c s ) is close to the quantum resistance [ItJ and therefore 
weakly depends on temperature. In our consideration we 
assumed, based on the experiments [H 0, 0, EH, that 
in the SFM state the sample resistance is much larger 
than the quantum resistance (l7j . meaning that it ex- 
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hibits variable range hopping behavior and therefore is 
more sensitive to temperature than the resistance of man- 
ganite systems. 

Recently the nanoscale granularity in manganese 
oxides was directly observed experimentally in 
La2-2a;Sri + 2KMn207, [29j]. The cluster structure in 
these perovskite materials is introduced by dopants, 
creating the individual weakly coupled nanodomains. 
To describe the MR in these materials one has to take 
into account electron localization within each cluster. 
This means that besides the tunneling conductance, gt, 
a finite grain conductance go has to be considered as 
well. In this case the total conductance can be written 
in the form g(T,m 2 ) = g (T) g t (m 2 )/[g (T) + g t (m 2 )}. 
Below the Curie temperature T c s , the tunneling con- 
ductance is small, <7t(l) -C go (T < T£), such that 
g(T < T c s , 1) = gt(l), whereas above T", the grain con- 
ductance go(T) becomes small due to localization effects 
(e.g. Jahn- Teller effect [1(3]), leading to the formation of 
an insulation state, see e.g. [H]. In our paper we were 
assuming that go ^> gt therefore the localization effects 
within each grain are small. This situation is realized in 
e.g. La^ArMnOa (A= Sr, Ca), P- 

The above considerations were carried out at zero ex- 
ternal magnetic field B. A finite field in ferromagnetic 
domain materials affects the SFM-SPM transition and 
leads to a reduction of the peak in the MR accompa- 
nied by a shift to higher temperatures with increasing B, 
which is parallel to the magnetization of the SFM at low 
T [Fig. Q] (bottom); Refs. [3, 0]. This means that the 
distribution function f(9) favors small angles near T c s , 
i.e. the drop in m 2 is smeared out. 

To summarize, we have investigated transport prop- 
erties of ferromagnetic nanoparticle arrays and nan- 
odomain materials in the limit of weak coupling between 
grains in the SPM and SFM states in both Ohmic and 
non-Ohmic regimes. We have described the electron 
transport near the Curie temperature in the artifi- 
cally self-assembled superspin arrays and discussed pos- 
sible applications of our results to a certain class of doped 
manganites, where localization effects within the clusters 
can be neglected. We derived the magnitude of the jump 
in the resistivity at the transition between the SPM and 
SFM states. We also discussed the influence of the mag- 
netic field on the jump amplitude and the relation of our 
results to available experimental data. 
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